Probable absence of a quadrupolar spin-nematic phase 
in the bilinear-biquadratic spin-1 chain 
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We study numerically the ground-state phase diagram of the bilinear-biquadratic spin-1 chain near 
the ferromagnetic instability point, where the existence of a gapped or gapless nondimerized quan- 
tum nematic phase has been suggested. Our results, obtained by a highly accurate density-matrix 
renormalization-group (DMRG) calculation are consistent with the view that the order parameter 
characterizing the dimer phase vanishes only at the point where the system becomes ferromagnetic, 
although the existence of a gapped or gapless nondimerized phase in a very narrow parameter range 
between the ferromagnetic and the dimerized regimes cannot be ruled out. 

PACS numbers: 75.10.Jm 
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I. INTRODUCTION 

In the past two decades, a large number of papers were 
devoted to the study of various properties of quantum 
spin chains. This was inspired to a large extent by Hal- 
dane's conjecture^ which states that isotropic antiferro- 
magnetic spin chains with half-odd-integer or integer spin 
values behave completely differently. While the excita- 
tions spectrum is gapless in the first case, a gap, the 
so-called Haldane gap, is generated in the other case. 

The most general form of isotropic coupling for S = 
1/2 spins is the usual Heisenberg model. For higher 
spin values, higher powers of the spins may appear in 
the Hamiltonian and this gives rise to a richer phase 
diagram. For 5 = 1 spins, assuming the most general 
isotropic nearest neighbor interaction, the Hamiltonian 
has bilinear and biquadratic terms in the spin operators 
and it can be written in the form 



i 

= [cos9(S t ■ S i+1 ) + sinfl(S 4 • S t+1 ) 2 ] 



(1) 



In this parametrization, 9 = corresponds to the an- 
tiferromagnetic model where the spectrum has a finite 
Haldane gap, while at 9 = ±7r the system is ferromag- 
netic. In fact both phases have a finite extension in the 
parameter space. The massive Haldane phase is stable 
for — 7r/4 < 9 < 7r/4, while ferromagnetism exists for 
7r/2 < 9 < it and —it < 9 < -37r/4. The phase bound- 
aries 9 = ±7r/4, — 37r/4 and 7r/2 are in fact special points, 
where the model can be solved exactly^S*^ 

The Haldane gap vanishes at 9 = it/ 4 and three soft 
modes appear at k = and ±27r/3. The system remains 
critical for it/ 4 < 9 < it/ 2 with the soft modes remaining 
at k = 0,±27r/3i^ii This gives rise to a power-law de- 
cay of correlations with a 3a periodic oscillation, hence 
the name "trimerized" phase, although the translation 
symmetry is not broken. 

The gap vanishes also at the other end of the Haldane 
phase, at — it/ 4, but it reopens for 9 < —tt/4 and a mas- 



sive dimerized phase with spontaneously broken transla- 
tional symmetry is found. A possible definition of the 
dimer order parameter D is 



D = lim 

N^oc 



\D 



N/2\ 



M+i/ 



(2) 



where chains are considered with open boundary condi- 
tion. Dimerization is measured as the alternation in the 
bond energy in the middle of a long enough open chain. 

The properties of this phase are best known for 9 = 
—it/2, where a partial mapping 8 to the 9-state quan- 
tum Potts model allows to calculate exactly the ground 
state energy, the gap, the correlation length, and also the 
dimer order^i 

While the properties and boundaries of the ferromag- 
netic, Haldane, and critical "trimerized" phases have 
been well established, there has been a long debate in 
the literature as for the boundary of the dimerized phase 
and the eventual existence of another phase between 
this dimerized phase and the ferromagnetic one near 
9 = —Stt/4, as seen in Fig. ^ 
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FIG. 1: Schematic plot of the phase diagram of the bilinear- 
biquadratic S = 1 model as a function of 6. 

By studying fluctuation effects near the end point of 
the ferromagnetic phase, Chubuko\ii£ claimed that there 
should be a gapped nondimerized nematic phase between 
the ferromagnetic and the dimerized phases. In fact, 
both ferromagnetism and dimerization involves sponta- 
neous symmetry breaking: for ferromagnetism SU(2) is 
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broken, whereas for dimerization it is translation invari- 
ance. These symmetries are largely unrelated, and there 
seems to be no a priori reason why the two should get 
broken hand in hand in one single transition. 

According to Chubukov's scenario the dimer order pa- 
rameter is finite in the dimerized phase, vanishes together 
with the gap at a 9 C close to, but definitely above — 3ir/4. 
The gap reopens for 9 < 9 C and closes again at — 37r/4, 
but the dimer order parameter remains zero in this whole 
range of 9. In this extended region, the system would 
have a nondegenerate singlet ground state and unbroken 
translational symmetry. Since the higher dimensional 
counterpart of Chubukov's phase would have quadrupo- 
lar order, this phase is usually called the "spin nematic" 
phase. The behavior of the gap and the dimer order pa- 
rameter according to this scenario is shown schematically 
in Fig. Eta). Other field theoretic calculations based on 
a nonlinear er-model approach for the director field also 
supported this scenario^ Note, however, that all these 
field theoretic calculations assumed a translation invari- 
ant ground state and did not check its stability against a 
possible spontaneous breaking of translation invariance. 
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FIG. 2: Schematic plot of the excitation gap (A) and the 
dimer order parameter D as a function of 9 around the phase 
boundary of the ferromagnetic and dimerized phases, (a) 
Chubukov's suggestion; (b) the result of earlier numerical 
work; (c) the case of a critical nondimerized phase. 

In our earlier numerical calculations?^^ where the 
vanishing of the gap was studied, we found no evidence 
for a closing and reopening gap. In these works both 
the gap and the dimer order parameter, indicating bro- 
ken translation invariance, were reported to vanish at 
9 C = — 37r/4 only. This scenario, the absence of the ne- 
matic phase is shown in Fig. Elk). 



Recent quantum Monte Carloi& and density-matrix 
renormalization-group (DMRG) 17 calculations have in- 
dicated that although Chubukov's proposal may not be 
completely correct in a closing and reopening gap, an 
exotic, critical phase with quadrupolar correlations may 
exist between the ferromagnetic and dimerized phases for 
— 37r/4 < 9 < 9 c . The behavior of the gap and the dimer 
order in this scenario is shown in Fig. Etc). Lauchli et 
alH estimated the value of 9 C to be — 0.677T. 

The aim of the present paper is to find further numer- 
ical evidence for the possible existence of this nondimer- 
ized phase. For this purpose we have studied finite spin 
chains up to 1000 lattice sites with open (OBC) and pe- 
riodic boundary conditions (PBC) for various 9 values 
using the DMRG method. 18 We have analyzed the behav- 
ior of the excitation gap and the dimer order parameter 
given by Eq. 



II. NUMERICAL PROCEDURE 

The numerical calculations were performed using the 
DMRG algorithm. Since the numerical accuracy is of 
crucial importance in the present study, this section is 
devoted to the problem of how the accuracy of our cal- 
culations could be determined and controlled. 

We have performed DMRG calculations both by us- 
ing the standard technique, 18 i.e., by keeping the num- 
ber of block states fixed and by using the dynamic block 
state selection (DBSS) approach ] 19 i 20 All eigenstates of 
the model have been targeted independently using two 
or three DMRG sweeps. 

In the standard procedure M = 500—1000 block states 
have been used. It was found that for the largest systems 
built up of N = 500 — 1000 lattice sites the truncation 
error varied in the range 10~ 8 — 10~ 9 for OBC and 10 -5 — 
10~ 7 for PBC. The following numbers are indicative of 
the accuracy: for OBC using M = 300 or M — 500 block 
states the ground state energies at 9 = —0.7ir differ in 
the fifth digit, ££(300,500) < 10" 5 , and the accuracy 
improved one order of magnitude when M = 1000 block 
states were kept, 5E(500, 1000) < 10~ 6 . 

The DBSS approach- 9 - allows for a more rigorous con- 
trol of numerical accuracy, and we set the threshold value 
of the quantum information loss \ to 10~ 8 . The mini- 
mum number of block states M m i n has been set to 256. 
The entropy sum rule was checked for all finite chain 
lengths for each DMRG sweeps, and it was found that 
the sum rule was satisfied after the second sweep already. 
The maximum number of block states varied in the range 
600 - 1400 for OBC and 1000 - 2500 for PBC, respec- 
tively. 

After accomplishing the infinite lattice procedure and 
using White's wave-function transformation methodSi 
the largest value of the fidelity error of the starting vector 
*stv, t>e* Btv = 1 - (*T|*stv) was of the order of 10~ 10 , 
where $t is the target state determined by the diagonal- 
ization of the superblock Hamiltonian. 
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As another test of the accuracy we calculated the dimer 
order profile Di using PBC. In theory, this should vanish 
identically for all finite chain lengths. Using the DBSS 
approach with \ = 10~ 6 , M m - m — 256 for chains up to 
N = 200 sites, in the parameter range — 0.757T < 9 < 
— 0.57T, the value obtained for Di was less than 10 -5 for 
all i = 1, . . . , N. The ratio of the DMRG block energy 
and the number of bonds within the block agreed up 
to 5 digits with the bond energy obtained between the 
two DMRG blocks. These results indicate that the finite 
dimer order parameter calculated with OBC is probably 
accurate at least up to 4 digits. 

III. GAP VS DIMER ORDER PARAMETER 

In order to obtain the energy gap A and the dimer 
order parameter D in the thermodynamic limit N — > oo, 
a finite-size scaling analysis has to be performed. In 
this section, as a benchmark case, we study in detail 
the exactly solvable point 9 = —tt/2. We demonstrate 
that for OBC, which is usually preferred in DMRG, the 
dimer order parameter is a better quantity to be ana- 
lyzed, as it provides much more accurate results. At the 
special point 9 = —tt/2 most of the quantities of spe- 
cial interest have been determined exactlyi&^iifliil The 
gap is A cxact = 0.173178, the dimer order parameter 
reads D exact = 1.124378, and the coherence length is 
£exact = 21.072. 
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FIG. 3: The excitation gap (A) and dimer order parameter 
(D) at 9 = -0.5tt as a function of 1/N 2 . The symbol x 
denotes the exact value. The solid lines are least square fits 
using the scaling forms discussed in the text. 

We used OBC and the DBSS approach with \ = !0~ 8 , 
Mmm = 256. Our results for the gap as a function of 1/N 
up to N — 500 are shown in the first panel of Fig. [3J A 
quadratic fit using the form in Eq. Q for 250 < N < 500 
yields A = 0.177 and a = 2300. The gap is about 2% 
higher than the exact result. 



A. Energy gap 

In a noncritical model with PBC the gap A(A^) is ex- 
pected to scale in leading order as 

A(A0=A + c-^expHV/£). (3) 

For OBC, however, the corrections are algebraic, and 
A(A^) is expected to vary as 

A(N) = A + a/N 2 + 0(N~ 4 ) (4) 

where a is a suitable constant. A qualitative argu- 
ment for this scaling ansatz can be given as follows^ 
The magnon dispersion is quadratic around its mini- 
mum, e(fc) = %/ A 2 + v 2 k 2 . Due to the boundary con- 
dition, however, the magnon wavefunction should have 
nodes on the boundary, which constraints the lowest 
possible magnon momentum to be k = ±tt/N. Conse- 
quently, the lowest possible excitation energy (the gap) is 
A(N) = y/A 2 +v 2 (tt/N) 2 m A + n 2 v 2 /2AN 2 + 0(N- 4 ), 
giving Eq. (@J. 

Since in the region of interest the ground state is a 
singlet and the lowest lying excited state is in the S to t = 2 
total spin sector, the excitation gap is calculated from the 
energy difference between the lowest lying levels of the 
Stot = 2 and Stot = sectors 

A(iv)=4°; t=2 (AT)-4°; t=0 (iv). © 



B. Dimer order 

For finite open chains with an even number of lattice 
sites the two typical valence-bond configurations of the 
dimerized phase are shown in Fig. 21 Due to the bound- 
ary condition these two singlet states are separated by a 
finite energy gap even in the thermodynamic limit. The 
ground state retains a high overlap with the configuration 
depicted in Fig. EJa) as N — > oo. This allows us to mea- 
sure the possible dimer order by considering D(N) and 
take the limit N — > oo. Note that this is only possible 
for OPB. For PBC the two possible dimerized configura- 
tions mix up and restore translation invariance for any 
finite N. In this case dimer order could only be measured 
by considering the dimer (4-point bond-bond) correlation 
function. In fact this was the quantity measured numer- 
ically in Ref. 0. Since calculations with PBC give less 
accurate results in DMRG than with OPB, we do not 
pursue here this approach. 

Accordingly, in our numerical calculations we used 
OBC with x - 10~ 8 , M min = 256 and determined D(N) 
as defined in Eq. J2J at the center of the chain. Two 
DMRG sweeps were taken and we checked that the en- 
tropy sum rule was satisfied. Our result is shown in the 
second panel of Fig. |3J The upward curvature of the data 
points as a function of 1/N is apparent for very short 
chain lengths (N ~ 40 - 80) and already for N > 200 
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FIG. 4: Typical valence-bond configurations for chains with 
even lattice sites. 

D(N) agrees up to three digits with the infinite chain- 
length limit. 

For noncritical models the characteristic behavior of 
the system is determined by a finite correlation length. 
Therefore, the end effects decay exponentially and the 
local quantity D(N) is expected to vary in leading order 
according to 

D(N) = D + dN- p e^>(-N/2^), (6) 

which is qualitatively similar to the PBC scaling of the 
gap in Eq. @, except that the scaling variable is the 
distance of the middle of the chain from the boundary, 
N/2, and the exponent of the algebraic prefactor is a pri- 
ori unknown. Nevertheless, knowing the exact value of 
D for 9 — — 0.57T, a least square fit provided /3 very close 
to 1. Using this, our numerical data for N > 60 can be 
fitted with D = 1.124375, £ = 20.2, and a = 5.9 (see the 
solid line in the second panel of Fig. |3J) . The correlation 
length is off by 4%, but the numerical value of D has an 
excellent relative accuracy of 3 x 10~ 6 . The exponential 
convergence of local quantities such as the dimer order 
parameter D(N) makes the extrapolation to the ther- 
modynamic limit very reliable. Our general conclusion 
is that — because of their different scaling behavior — the 
dimer order parameter is a much better quantity to ana- 
lyze than the energy gap. In the next section we pursue 
this idea to investigate the phase diagram as a function 
of 9. 



IV. NUMERICAL RESULTS 

For the reasons presented above, the highest chance 
to find the subtle nondimerized spin-nematic phase near 
9 = — 37r/4 is by studying the dimer order in open chains. 
Therefore, we calculated D(N) using OBC up to 1000 lat- 
tice sites for various 9 values between — 37r/4 and — n/2. 
We have set x — 10 -8 , M m i n = 256 and used three 
DMRG sweeps. Due to the very large correlation length 
in the vicinity of the ferromagnetic phase boundary the 
maximum value of M is varied in the range of 600 — 1400 
for our longest chains. Our numerical results are shown 
in Fig. |51 for two different 9 values. The numerical error 
is much smaller than the size of the symbols. 

It is seen in the figure that for large N the data points 
show an upward curvature as a function of 1/N for all 9 
values. On the other hand, the inflexion point shifts to- 
wards very long chain lengths with decreasing 9 values. 



9 = -0.6731 
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FIG. 5: Finite size scaling of the dimer-order parameter at 
9 = -0.67tt and -0.70tt for 100 < N < 1000. 



The extrapolated value D of the dimer order parameter 
was determined using Eq. 0. These values are shown 
in Fig. [S] as a function of 9. In fact, we found a finite, 
nonzero D for all 9 shown. However, we have no results 
for 9 < 0.77T, where D(N) is so small that it is comparable 
to the numerical error. Although the dimer order param- 
eter decreases very rapidly for 9 < — 0.64-7T, the smooth 
behavior of D as a function of 9 and the upward curva- 
ture observed suggests that it vanishes at 9 — — 3ir/4. 

It is also apparent from Fig.^that the dimer order pa- 
rameter resembles the form of the Berezinskii-Kosterlitz- 
Thouless (BKT) transition, thus it opens exponentially 
slowly as a function of 9, 

D(9) = aexp[-c(9-9 c )-°}, 9 c = -3tt/4, (7) 

for 9 > 9 C . In Eq. l|7|. a and c are nonuniversal constants 
and a is a characteristic exponent, which was estimated 
by a least square fit to be a = 1.3 ± 0.3 (see Fig. EJ. 
This functional form is in qualitative agreement with our 
earlier numerical result^* 

Based on our calculations we conclude that there is 
no sign of any phase transition to either a gapped or a 
gapless nondimerized phase in the vicinity of 9 ~ — 0.677T, 
and thus a direct phase transition takes place between the 
ferromagnetic and the dimerized phases. 



V. CONCLUSION 

In summary, we have performed a density-matrix 
renormalization-group calculation on the spin-1 bilinear- 
biquadratic spin chain model in the vicinity of the fer- 
romagnetic phase in order to search for a nondimerized 
quantum nematic phase suggested by Chubuko\ii£ or an 
extended critical region reported recently by Lauchli et 
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FIG. 6: The extrapolated value of the dimer order parameter 
as a function of 6/tt. The solid curve is a least square fit using 
the form in Eq. with a = 2.4 ± 0.6, c = 0.11 ± 0.05, and 
a = 1.3 ±0.3. 

aJJi We took special care of the numerical accuracy since 
it has special importance in the present problem. We 
used the DBSS approach with the maximum number of 
block states varying between 1000 and 2000, and per- 
formed calculations on very long chains up to N = 1000 
lattice sites. 

As a benchmark case we analyzed in detail the ex- 
actly solvable point 9 = —tt/2. We have found that the 
dimer order parameter is a quantity, which can be deter- 
mined numerically much more accurately than the energy 
gap, because of their unequal finite-size scaling forms. 
Whereas the gap scales algebraically in open chains, the 
dimer order parameter, which is a local quantity mea- 
sured in the middle of the chain, becomes highly inde- 
pendent of end effects and scales exponentially. 



The phase diagram of the model was explored in gen- 
eral by computing the dimer order parameter as a func- 
tion of 9. We have found strong indications that the 
dimer order, which characterizes the dimer phase, only 
vanishes at the phase boundary of the ferromagnetic 
phase 9 = —Sir/A. The behavior resembles that of 
Berezinskii-Kosterlitz-Thouless transition, i.e., the order 
parameter opens exponentially slowly as we move away 
from the transition point. Our findings are in agreement 
with the exponentially slow opening of the energy gap re- 
ported earlier M Nevertheless, the actual transition itself 
is of first order as it involves direct level crossings, as well. 
It is noteworthy that at 9 = — 37r/4 the model has an ex- 
tra SU(3) symmetry^ which causes extra degeneracies 
in the spectrum. We are tempted to speculate that this 
extra symmetry may play a role in that the two phases 
with seemingly unrelated broken symmetries, SU(2) and 
translation invariance, meet in this special point without 
an intermediate phase. 

We could not find any trace of a nondimerized regime, 
at least surely not above 9 ss 0.77T. Below this limit nu- 
merical precision is a crucial issue since the quantities of 
interests are extremely small. If the intermediate phase 
exists, it should be constrained in a very narrow region 
near 9 = — 37r/4, certainly much narrower than predicted 
by Lauchli et ai. The more likely alternative interpre- 
tation, i.e., the nonexistence of the intermediate phase, 
is clearly at odds with current field theory analysis. A 
reconciliation of the numerical results with field theory 
would be very welcome in the future. 
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